A quantum criticality perspective on the charging of narrow quantum-dot levels 
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Understanding the charging of exceptionally narrow levels in quantum dots in the presence of 
interactions remains a challenge within mesoscopic physics. We address this fundamental question 
in the generic model of a narrow level capacitively coupled to a broad one. Using bosonization 
we show that for arbitrary capacitive coupling charging can be described by an analogy to the 
magnetization in the anisotropic Kondo model, featuring a low-energy crossover scale that depends 
in a power-law fashion on the tunneling amplitude to the level. Explicit analytical expressions for 
the exponent are derived and confirmed by detailed numerical and functional renormalization-group 
calculations. 
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Introduction. Confined nanostructures offer a unique 
arena for thoroughly interrogating the interplay between 
interference and interactions while holding the promise of 
future applications. Particularly appealing are semicon- 
ductor quantum dots (QDs), for which the manipulation 
of spin [l|, [IJ and charge Q has recently been demon- 
strated. The precise and rapid control of switchable gate 
voltages renders these devices attractive candidates for 
a solid-state qubit [4, 5]. The accurate manipulation of 
QD setups requires, however, detailed understanding of 
how charging proceeds. Indeed, interactions can substan- 
tially modify the orthodox picture of charging, whether 
by renormalizing the tunneling rates or by introducing 
nonmonotonicities into the population of individual lev- 
els [(| 0, HI- Even the simplest two-level device, where 
each level harbors only a single spinless electron, displays 
remarkably rich behavior 

We consider a situation in which the width of one nar- 
row level is much smaller than the width of the other 
broad one. A disparity in widths is generic for QDs in the 
intermediate regime between integrable and chaotic Q. 
It was reported in several artificial structures 1(| U | , and 
has been exploited for charge sensing 

BUI- 

As the en- 
ergy e_ of the narrow level is raised, its occupation varies 
from 1 to over a characteristic width fi. This energy 
scale, or the corresponding charge-fluctuation time scale 
H/Cl, manifests itself in charge sensing and transmission- 
phase measurements [3] ■ The effect of inter-level repul- 
sion U on SI has been explored only in the larg e- U limit, 
revealing novel correlation effects [l5|, [l(| llTl llij]. The 
physical mechanism determining SI for moderate U re- 
mains unclear 0]. 

In this Letter we solve the fundamental question of the 
charging of a narrow QD level from a quantum-critical 
perspective. Due to the capacitative coupling U, every 
switching of the narrow level initiates restructuring of the 



broad level and its attached Fermi sea, in direct anal- 
ogy with the x-ray edge singularity. For nonzero tunnel- 
ing to the narrow level, coherent superpositions of these 
charge re-arrangements lead to Kondo physics (l9| with 
the charge state (0 or 1) acting as a pseudo-spin, and 
the energy of the narrow level acting as a Zeeman field. 
Using Abelian bosonization we show that SI, being the 
Kondo scale in the pseudo-spin language, depends on the 
tunneling amplitudes in a power-law fashion. We derive 
explicit analytical expressions for the exponents encom- 
passing all physical regimes of the model (at zero temper- 
ature T) . In a second step we confirm our predictions by 
detailed numerical ^ (NRG) and functional 2l| (FRG) 
renormalization-group (RG) calculations, thus resolving 
this challenging aspect of mesoscopic physics. 

Model and objective. Our specific model for charging 
is depicted schematically in the inset of Fig. 1, and is 
defined by the Hamiltonian (a is the pseudo-spin index) 
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Here, d± (cj. ± ) creates an electron on the dot (in the 
leads), and An± equals d±d± — 1/2. Equation |T]) is a 
generalized Anderson impurity model with pseudo-spin- 
dependent tunneling amplitudes V + > V_ > and a 
tilted magnetic field, whose components are e + — e_ and 
the direct hopping amplitude b. This form follows from 
a generic model of spinless electrons with two dot levels 
and two leads by simultaneous unitary transformations 
in the dot and the lead space 

[II ELI El. Th e Hamil- 
tonian (TT]) has recently gained considerable attention in 
connection with phase lapses, population inversion, and 
many-body resonances Q . The energies e ± are tuned us- 
ing gate voltages. Depending on the specific realization, 
their tuning may inflict a similar change in b. We focus 
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FIG. 1: (Color online) The exponent a computed using the 
NRG, FRG, and Eq. ©. NRG parameters: F+/D = 0.04, 
A = 1.7, and 2800 states are retained. Inset: The model 
system. Two localized QD levels are coupled by tunneling to 
separate baths. Spinless electrons residing on the two levels 
experience a Coulomb repulsion U. 



on realizations where e_|_ can be tuned independently of 
b. 

The bare energy scales that characterize tunneling in 
Eq. ([1]) are the level broadenings T± = irpV± and the 
direct hopping amplitude b. The density of states (DOS) 
p is taken to be equal for both bands without loss of 
generality. Our interest is in the charging properties of 
the narrow level d_ as a function of e_ in the limit where b 
and V_ are both small: r_, b < T + . Strictly at b=V_=0 
ergodicity of the microcanonical ensemble is broken as a 
new conserved quantity arises: n_=d_d_ is either equal 
to or 1. Comparing the total energies of the competing 
ground states with (n_) — and (h_) — 1 as a function 
of €_ one finds a critical value e_ = e*(e+, U, V+) at which 
the two become degenerate. For e + = 0, e* is pinned to 
zero by particle-hole symmetry if symmetric bands are 
assumed. In the limit b, V— — > the average occupation 
(n_) thus indicates a first-order transition (width fi = 0) 
as e_ is swept across e*. It is the smoothening (SI > 0) 
of this transition at small but finite b, V_ for U > that 
is addressed in this Letter. 

Two regimes can be distinguished depending on e + . 
When |e + | ^> U, T + , the level d} + maintains an approxi- 
mately fixed integer valence (n + ) G {0,1}, independent 
of e_. Hence, the charging of d_ is essentially single- 
particle in nature with il = T_ + T + b 2 /e 2 ^. The effect of 
interactions is contained in the simple Hartree renormal- 
ization, e_ — > e_ + U((n + ) — 1/2). Far more complex is 
the case of |e + | <C max{[/, when the broad level is 
prone to strong valence fluctuations (for e_ — > e*). Going 
from U/T + < 1 to 1 < U/T + spans all physical regimes 
from weak to strong electronic correlations \u[ [ll| , 
which constitutes the main focus of our study. To this 
end we initially set e + = 0, which fixes e* — 0. Using an- 
alytical and numerical tools we first obtain Q in the case 
where either V_ or b is nonzero. The combined effect of 
V_ and b is next addressed by single-parameter scaling 



and FRG. Finally, we extend our analytical results to 
arbitrary e+. 

Analytical approach. To analytically determine the 
width Q using minor approximations, we proceed in two 
steps. First, we derive a continuum-limit Hamiltonian 
where r + is incorporated in full. Second, an exact map- 
ping of this Hamiltonian onto the anisotropic Kondo 
model is established. This allows usage of known results 
for the Kondo problem in order to extract SI. 

In the first step, we diagonalize the Hamiltonian 

= T,kVl+ c k+ + v +T,k{ c l+ d + + d + c k+] usin S 
scattering theory. Expanding in terms of the single- 
particle eigen-modes of and converting to continuous 
constant-energy-shell operators [13] , H takes the form of 
a generalized interacting resonant-level model with a sin- 
gle d_ level tunnel coupled to two bands: a narrow a = + 
band with a Lorentzian DOS of half- width T + , and a flat 
a = — band with half- width D 3> T + . In addition, the 
d_ level is capacitively coupled to the '+' band. 

In the desired limit b, T_ <C T + , one can conveniently 
replace the Lorentzian DOS with a flat symmetric one 
of height l/7rr + and half- width D + = nT + /2 22]. The 
elimination of all degrees of freedom in the energy interval 
D + < |e| < D leads to renormalizations of the couplings 
of the order of r_/T + < 1 or higher, which can be safely 
neglected. Converting at this point to left-moving fields, 
we obtain the continuum-limit Hamiltonian 

/oo 
, ipl.(x)d x 'ip a (x)dx + e_d_d_ 
-oo 

+ (6/2)Va{V , + (°) d - + H - c -} + Ua :V'+(°)^+(°) :A ^- 
+ v /ar+r_{?/'i(0)d_ + H.c.} , (2) 

applicable at energies below T + . Here, a = ■kTivf/D + 
is a new short-distance cutoff ("lattice spacing"), and 
: ip\_ip+ : stands for normal ordering with respect to the 
filled Fermi sea. The left-moving fields obey canonical 
anticommutation relations subject to the regularization 
(5(0) = 1/a. The derivation of Eq. @ is controlled by the 
small parameters r_/T + <C 1 and b/T + <C 1, and hence 
is expected to become asymptotically exact as T_ , b — > 0. 

If either b = or T_ =0, Eq. ^ can be treated using 
Abelian bosonization [23| . To this end, we introduce two 
bosonic fields $>±(x), one for each fermion field ip±{ x )- 
With a proper choice of the phase-factor operators, the 
bosonized Hamiltonian reads 

f- /-OO 

H = Y — [V$ (T (x)fdx + e_dld_ 
+ hv F ^V®+(0)An_ + A/ e l *±(°)d_ + H.c.j, (3) 

7T V2 L J 

The tunneling term in Eq. ([3]), proportional to A, de- 
pends on the case of interest; one takes A = 6/2 and the 
upper sign (A = v/T^TT, lower sign) for r_ = (6 = 0). 
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The value of Sjj = arctan(£//2r + ) is fixed by matching 
the b = T_ = scattering phase shifts of the '+' band in 
the fermionic and the bosonic representations, for each 
sector with fixed integer occupancy of the '— ' level. 

Next, we manipulate Eq. ([3]) by (i) applying the canon- 
ical transformation TV = WHU with 



U = exp [-«(2<St//7r)$+(0)An_] 



(4) 



and (ii) converting to the "spin" and "charge" fields 
$> s (x) and Q c (x)- The latter are defined as & s (x) = 
$+(x) and $ c (x) = $_(x) for T_ = 0, and 
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$ T (x) T — *±(a0 (5) 
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for 6 = (the upper signs correspond to $ s ). In this 
manner, the Hamiltonian acquires the unified form 
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where 7 = + (25u/t) 2 for 6 = and 7 = 1 — 2Su/tt 
for r_ = 0. 

The very same Hamiltonian with < 7 < \[2 also 
describes the anisotropic Kondo model with < J z , 
where in standard notation A = J±/V& and 7 = \/2~[l — 
(2/tt) arctan(7rpJ z /4)] represent the transverse and lon- 
gitudinal spin-exchange couplings respectively, and e_ = 
(iBgH corresponds to a local magnetic field. This rep- 
resentation of the Kondo model is obtained by [24| (i) 
bosonizing the Kondo Hamiltonian with two bosonic 
fields $t( :E ) an d < J ) l( a; ); (ii) converting to the spin and 
charge fields $ SlC (x) = [<&i(x)=F$i(x)]/\/2, (iii) employ- 
ing W = fiHT with f = exp[-i\/2(2* je /7r)$ 8 (0)r 2 ], r 2 
being the z spin component and 5 Z = arctan(7T/9j z /4), 
and (iv) representing the spin r in terms of the fermion 
d_ = r~ . This establishes a mapping between our prob- 
lem with either b = or T_ = and the anisotropic 
Kondo model. In particular, charging of the d_ level is 
mapped onto the magnetization of the Kondo impurity, 
relating the width SI to the Kondo temperature 2V • 

We can now exploit known results for the Kondo 
problem. Specifically, RG equations perturbative 
in J± but nonperturbative in J z [T3] give Tr ~ 
D + (A/D + ) 2 /<- 2 -> >, which yields for our problem 
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(b/T+fP if r_ = , 
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Thus, is a power law of the relevant tunneling ampli- 
tude with an exponent that varies smoothly with U. In 
going from U = to U 3> T + , a grows monotonically 




FIG. 2: (Color online) The exponent /3 computed using the 
NRG, FRG, and Eq. ©. NRG parameters: F+/D = 0.02, 
A = 1.6, and 2000 states retained. Inset: representative NRG 
data for Q vs b, along with the log-log fits used to extract j3. 



from 1 to 7rf//(8r + ) while (3 decreases from 1 to 1/2. 
The asymptote a — irU/(8T + ) coincides with the result 
of Ref. 17 [Eq. (29) with e = -U/2], obtained using 
very different techniques. For U = 0, the noninteracting 
integer exponents are reproduced. Hence Eqs. (JSj) are 
precise both at small and large U. As shown next, these 
expressions remain highly accurate also at intermediate 
U, suggesting that they might actually be exact. 

Numerical analysis. To test Eqs. ([5J, we computed a 
and j3 numerically using the NRG [20( and FRG 
each approach having its own distinct advantage. 
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NRG is extremely accurate in all parameter regimes of 
interest, while the FRG is approximative in U but offers 
a far more flexible framework for scanning parameters. 
The width Q = l/ (kXc) w as obtained with either method 
from the inverse charge susceptibility Xc — d(n_) /de_, 
evaluated at e_ = and T — > 0. The exponents a and /3 
were extracted from log-log fits (see the inset to Fig. [2]) . 
Our results, summarized in Figs. [TJ and [51 reveal excel- 
lent agreement between Eqs. © and the NRG, to within 
numerical precision. The agreement extends to all inter- 
action strengths from small to large U, confirming the 
accuracy of Eqs. ([8|) at all U. The FRG results for a 
coincide with those of the NRG up to U/T + w 2, above 
which they acquire a linear slope that is reduced by a 
factor of 8/71" 2 as compared to the NRG (l7j |. The ex- 
ponent (3 is accurately reproduced up to larger values of 
U/T + . In particular, the FRG data for a and (3 exactly 
reproduce the leading behaviors of Eqs. ((8]) at small U. 

Combination o/T_ and b. The case where both F_ and 
b are nonzero lies beyond the scope of our bosonization 
treatment, but allows the formulation of a scaling law. To 
this end, consider the dimensionless quantity SI = Q/ D + , 
which depends on the three dimensionless parameters in 
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f(y,b,6 v ), with v = vTvrr/^4 



and 



b = b/D + . Given the exact RG trajectories, evolves 
according to ST = fi/f = f(V', b', S' v , {A^}) upon re- 
ducing the bandwidth from D+ to £D + (0 < £ < 1). 
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FIG. 3: (Color online) A scaling plot of fl for two (fixed) ra- 
tios U/T+ and different combinations of T_ and b, obtained 
using the FRG. The coefficients A and B were extracted from 
the limiting cases where 6 = and F_ = 0, respectively [25| . 
For clarity, the data for U /r+ = 5 were multiplied by a con- 
stant as indicated by the arrow length. Dashed lines show the 
asymptotes T = 1 and T — x. 

Here, primes denote renormalized parameters and {A^} 
are the new couplings generated. At sufficiently weak 
tunneling the RG equations can be linearized with re- 
spect to the relevant couplings V' and I >', resulting in 
their power-law growth with the exponents determined 
previously: V' = V£-V 2a and V = b^ 1/2p . Note 
that 5u is left unchanged in this approximation, nor 
are there any new couplings generated. Consequently, 
f(V, b, Su) = £f(Vr 1/2a , ir 1/2p , Su) = Q is a homoge- 
neous function of £, taking the general form f{V, b, Sjj) — 
V 2a g{b 2 P /V 2a ,Su). Finally, defining the coefficients A 
and B from f2|&=o = AT™ and S7|r_ =n = Bb 213 , we arrive 
at the scaling form (25j 

n = AT™F(Bb 2 < 3 /AT™;5 u ), (9) 

with F(0;6u) = 1 and T(x > l;Su) = x. In Fig. [3] we 
confirm the scaling form of Eq. ^ using FRG data. 

Extension to arbitrary e. . Our discussion has fo- 
cused thus far on e + = 0. A nonzero e + introduces the 
potential-scattering term H ps = e + a : ip^ + (0)ip+(0): into 
Eq. @. Consequently, Su in Eq. ^ is replaced with 
two distinct parameters S± = arctan[(C7 ±2e + )/2r + ], as- 
signed to An_ = ±1/2, respectively. An identical deriva- 
tion, only with 25jjh_ — > (S + + i5_)An_ + (5+ — 6-)/2 in 
Eq. ([4]) , leads then to the same Hamiltonian ([6]) with two 
modifications: (i) e_, and thus e*, acquires a shift propor- 
tional to 5+ — S 2 ^, and (ii) b~u is replaced with (S + + S-)/2 
in the expressions for 7. The end results for a and (3 
are just Eqs. © with 5u — * (S+ + S-)/2, which prop- 
erly reduce to the noninteracting limit a = (3 = 1 when 
|e + | ^> U,T + . The effect of nonzero e + is negligible for 
|e, I <C max{[/, r + }. It becomes significant only as |e, | 
approaches max{?7, 

Summary. We have resolved the fundamental question 
of the charging of a narrow QD level capacitively coupled 



to a broad one. The zero-tunneling fixed point is criti- 
cal in the sense of being unstable. Finite tunneling is 
a relevant perturbation, driving the system to a strong- 
coupling Fermi- liquid fixed point. The inverse charge- 
fluctuation time f2 varies as a power of the bare tunneling 
amplitude, with a nonuniversal exponent that depends 
on the nature of tunneling, the strength of the capaci- 
tive coupling, and the width and position of the broad 
level. We have proven this scenario by devising a two- 
stage mapping of the original model onto the anisotropic 
Kondo problem, yielding accurate analytic expressions 
for the exponents. Our analytic predictions were con- 
firmed by extensive numerical calculations within the 
frameworks of the NRG and FRG. 
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